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Phase-equivalent transformation of local interaction is generalized to the multi-channel 
case. Generally, the transformation does not change the number of the bound states in the 
system and their energies. However, with a special choice of the parameters, the transforma- 
. . . tion removes one of the bound states and is equivalent to the multi-channel supersymmetry 

^^ I transformation recently suggested in |g^. Using the transformation, it is also possible to 

^D . add a bound state to the discrete spectrum of the system at a given energy E < ii the 

^—^ ' angular momentum at least in one of the coupled channels I > 2. 

(N ■ 

<U ■ I. INTRODUCTION 

^: 

^ ■ Nucleon-nucleon, nucleon-cluster and cluster-cluster potentials are an input for various microscopic 
i calculations of nuclear structure and reactions. Unfortunately, the exact form of the potentials 
'""' [ describing these interactions is unknown. It is conventionally supposed that the interactions are 
^ ■ local, that is, of course, an approximation only. However, the available scattering data and bound 
0\ i states properties can be fitted with approximately the same accuracy by different local potentials. For 
example, there is a lot of so-called realistic A^A^ potentials on the market describing A^A^ scattering 
^' ■ and deuteron properties with high accuracy. More, a description of phenomenological data can be 
O ■ achieved with the potentials very different in structure. In particular, meson-exchange A^A^ potentials 
like the Nijmegen one |]1[, are known to have a short-range repulsive core in a triplet s wave. The 
same high-quality description of the nucleon-nucleon data is provided by latest versions of Moscow 
potential 0,^] that does not have a repulsive core but instead is deeply-attractive in the triplet s 
wave at short distances and supports an additional forbidden state. The possibility of alternative 
p] . description of various cluster-cluster and nucleon-cluster interactions by means of repulsive-core and 
deeply-attractive potentials with forbidden states, is also well-known (see, e.g., the discussion in p| 
and references therein). 

Principally it is possible to distinguish experimentally between alternative potentials studying their 
off-shell properties in interaction with an additional particle. The simplest probe is the photon, and 
as it was shown in 0-0, the proton-proton bremsstrahlung reaction pp —>■ pp^ at the energy range 
of 350-400 MeV can be used to discriminate between various nucleon-nucleon potentials. However 
the pp — »• pp^ reaction has not been examined experimentally in this energy range. 

Another possibility is to study properties of three and four body systems bound by two-body 
potentials of interest. From this point of view, it looks like that we do not have at present satisfactory 
nucleon-nucleon, cluster-nucleon and cluster-cluster potentials. It is well-known that none of the 
realistic A^A^ potentials provides proper binding of tritium or '^He. There are successful attempts in 
generating phenomenological three-nucleon interactions tuned to fit the properties of light nuclei 
(see also ^ and references therein). However, as it was shown in a detailed study of Picklesimer et 
al [^ , the effect of three-nucleon forces consistent with realistic two-body ones on the binding energy 
of the triton is canceled by effects of the virtual excitation of A isobars, etc. Hence the trinucleon 
cannot be satisfactorily described using known realistic two-body potentials supplemented by three- 
body potentials consistent with them. All calculations within three-body cluster models also fail to 



O 
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reproduce the correct binding energy of three-cluster nuclear systems with known local cluster-cluster 
and cluster-nucleon potentials fitted to the corresponding scattering data. 

To design a potential consistent with two-body phenomenological data and providing the correct 
binding of few-body systems, it seems promising to make use of phase-equivalent transformations 
depending on a continuous parameter (s). Some attempts in this direction have been performed using 
non-local phase-equivalent transformations. The results of these attempts are encouraging: in Ref. 
|T0| an oversimplified A^A^ potential providing a satisfactory description of s wave A^A^ scattering date 
was fitted to reproduce exactly the triton binding energy, while in Ref. [|1T| realistic n-a potentials 
were tuned to reproduce various ^He properties including the binding energy within the a + n + n 
cluster model. The interactions suggested in Ref. |]TU|,|n| are non-local. Various applications (see, 
e.g. 
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of local phase-equivalent transformations to few-body problems were restricted to the 
supersymmetry transformation P^^TB] that removes one of the bound states in a two body system. 
The supersymmetry transformation does not contain parameters and cannot be used for fine tuning 
of the interaction of interest. 

A local phase-equivalent transformation which preserves the number of the bound states and 
depends on a continuous parameter, exists and is well-known in the inverse scattering theory [|17 



studied in Ref. fl8 



Recently the effect of this transformation on the properties of three and four nucleon systems was 

using as an example a semi-realistic Malfliet-Tjon A^A^ interaction [0. It 
that a slight phase-equivalent modification of A^A^ interaction is enough 



was shown in Ref. [|T^ 



to reproduce the trinucleon binding energy and to improve simultaneously the description of four- 
nucleon binding. However the local transformation was developed for a single-channel case only and 
cannot be applied without some approximations to realistic A^A^ interactions that mix triplet s and 
d partial waves. Another drawback of the transformation is that it involves a bound state wave 
function and thus cannot be used to modify nn and pp interactions and the np interaction in all 
'non-deuteron' partial waves. 

Recently Sparenberg and Baye |2y] suggested a multi-channel supersymmetry transformation. We 
use some ideas of Ref. |^0|] to derive in what follows a multi-channel phase-equivalent transformation 
which depends on continuous parameters. The transformation can be treated as a generalization 
both of the single- channel phase-equivalent transformation |T^ and of the multi-channel supersym- 
metry transformation of Ref. [^. Generally, the transformation does not change the number of the 
bound states in the system and their energies. However, with a special choice of the parameters, 
the transformation removes one of the bound states and becomes equivalent to the multi-channel 
supersymmetry transformation suggested in ||20|. If the angular momenta in all coupled channels are 
less than 2, a parameter-dependent family of local interactions phase-equivalent to the given initial 
one can be constructed by means of the transformation even in the case when the system does not 
have a bound state. If the angular momentum at least in one of the coupled channels / > 2, the 
transformation can be used to add a bound state to the discrete spectrum of the system at a given 
energy E < 0. Having a bound state, one can construct a family of phase-equivalent potentials and 
afterwards remove the bound state by the supersymmetry version of the transformation. Thus, the 
suggested transformation can be used in a multi-channel case to produce phase-equivalent interac- 
tions without any restriction on the structure of the discrete spectrum of the system. In particular, 
the transformation can be applied to the realistic A^A^ interaction in all partial waves. 

II. GENERAL FORM OF LOCAL MULTI-CHANNEL PHASE-EQUIVALENT TRANSFORMATION 



Multi-channel scattering and bound states we describe by Schrodinger equation 



J2{H^,-E6,,)^,{E,r) = 0, 



where indexes i and j label channels, E is the energy, the Hamiltonian 
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dr"^ r 
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S^J + V,,{r), (2) 



m is the reduced mass, and k stands for the angular momentum in the channel i. We suppose that 
the potential Vij{r) (i) is Hermitian and (ii) at large distances it tends asymptotically to a diagonal 
constant matrix, 

Vij{r) — *€i6ij, (3) 

r— >oo 

where e^ is a threshold energy in the channel i. We suppose that ei = and e^ > ej ii i > j. 
Boundary conditions for the wave functions are 

ip,{E,0) = 0, (4) 

ifii^E, oo) < oo. (5) 

Except for the discussion in section piQ we suppose that there is at least one bound state in the 
system at the energy Eq. The corresponding wave function, (fii{Eo,r), is supposed to be normalized. 



oo 

^ / ifi*{Eo,s) 'fi{Eo,s) ds = 1, 



(6) 



where * denotes the complex conjugation. Of course, Lpi{EQ,r) fits more severe boundary condition 
at r — i> oo than (j^): 

(/^,(£;o,oo) = 0. (7) 

We define the transformed potential Vij{r) as 

V^J{r) = V.J{r) + v.,ir), (8) 

where 

„,,(.) ^-2cf^ ^.(Eo.r)^/(g„.r) _ ^^^ 

k a 

and y4, C and a are arbitrary real parameters. 
The main result of this paper can be formulated as the following statement. 

• The wave function 

r 

Y, f ipk*{Eo, s) (pk{E, s) ds 

^,{E,r) = ^,{E,r) - Cv,{Eo,r)±^ (10) 

A + CZl\MEo,s)\' ds 

k a 

Gts inhomogeneous multi-channel Schrodinger equation 



n^ 



J2{H,,-ES,,)^,iE,r) = C — 



'^i{Eo,r) 



where the Hamiltonian 



H^3 Srj 2^ 



h^ r d^ ^kih + iy 



dr"^ r^ 



V^A^) 



W{Eo,E;a) 



A + Cj:j\MEo,s)f ds 

k a 



and the quasi-Wronskian 

W{Eo, E;a) = Y, [Vk*{Eo, a) ^k'{E, a) - ^k*'{Eo, a) ME, a)] . 



;ii) 



(12) 



(13) 



We use prime to denote derivatives: /' = -^f. 
To prove the statement, one can verify Eq. (^) by the direct calculation of 

J2j [Hij — ESij] Lpj{E,r) using the definitions (|8|)-(|lOD, (|T^ and (|13|) and other formulas given 

above as well as the fact that the interaction Vij{r) is Hermitian, V*Ar) = Vji{r). The calculation is 
lengthy but straightforward. 

It is clear from (|1^) and (0) that the suggested transformation is phase-equivalent at any energy 
£" > 0; all the bound states supported by the initial potential Vij are preserved by the transformation 
since the wave functions (pi{Eb,r) for the corresponding energies Eb < (including Eq) fit both 
boundary conditions (^) and (|^). However, the denominator in the last term in (|I0|) should be 
non-zero at any distance r, and therefore one should be accurate in assigning values to arbitrary 
parameters A, C and a. This requirement can be easily satisfied in a wide and continuous range of 
parameter values. 



III. PARTICULAR CASES OF THE PHASE-EQUIVALENT TRANSFORMATION 
A. Homogeneous Schrodinger equation 

Of course, we are mostly interested in phase-equivalent transformations that results in homoge- 
neous Schrodinger equation 



Y,{Hij-E5i^^j{E,r) = Q 



(14) 



instead of the inhomogeneous Schrodinger equation (|ri|). To derive the transformation leading to 
Eq. (p!^, we can fix the parameters A, C, and a in such a way that the r.h.s. of Eq. ( [Tl| ) will take zero 
value. The choice C = brings us to the equivalent (contrary to phase-equivalent) transformation 
that is of no interest. Thus we should search for the parameters that fit the equation 



W{Eo,E-a)=Q. 



(15) 



Two obvious solutions of Eq. (|I5|) are a = and a = oo. Various other solutions of Eq. (|1^) can be 
found for particular potentials Vij{r). However, the non-zero finite solutions a of Eq. (|r^) are energy- 
dependent. With the solutions a{E) of Eq. (|1^) we can obtain energy-dependent potentials Vij{E] r) 
phase-equivalent to the initial energy- independent potential Vij{r). It may be interesting for some 



applications, but we shall not discuss the energy-dependent transformation and shall concentrate 
our attention on the solutions a = and a = oo. 
The case a = presents a generalization of the single-channel phase-equivalent transformation 



of Ref. 1 17]. For the bound state at the energy Eq, the wave function obtained by means of the 



transformation is of the form: 



iPi{Eo, r) = . (16) 

A + Cj:j\cpj{Eo,s)f ds 
j 



The wave function (|T^) is not normalized. The normalization constant can be easily calculated. The 



normalized bound state wave function is 



A + C..J, . VA{A + C)^,{Eo,r) 

-—— ipi{Eo, r) = . (17) 

A + Cj:f\^,{Eo,s)f ds 
j 

It is interesting that the components of the bound state wave function in all channels are modi- 
fied by the transformation sinchronically: all the components ipi{EQ,r) are multiplied by the same 

multiplier ^/A(A+~C)i A + C^ J \(pj{Eo, s)\ ds\ . Nevertheless the relative weight of the com- 

V i / 

ponents ifi{Eo,r) in the norm of the total multi-channel wave function can be changed by the 
transformation. 
Now let us discuss the case a = oo. The transformed wave function in this case is of the form: 

r 

Cipi{Eo, r) ^ / ipj*{Eo, s) ipj{E, s) ds 

UE, r) = ^.{E, r) '-^ . (18) 

A + CY.j\vAEo.s)\' ds 

j oo 

\i E ^ Eq, the functions i^i^E, r) and (pi{Eo, r) are orthogonal: 

ipC{EQ,s)ipi{E,s)ds = Q. (19) 





With the help of (|1^) and (^, we can rewrite (0) as 



Cipi{Eo,r)Y.j<^j*{EQ,s) iPjiE,s) ds 
i_o 

A-C + Cj:j\^,{Eo,s)f ds 



^,{E, r) = if,{E, r) '-^—^ . (20) 
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It is seen from ( pO]) that the case a = oo is identical (up to the redefinition of the parameter 
A ^ A + C) to the case a = ii E ^ Eq. It is clear, however, that after the redefinition of the 
parameter A ^ A + C, the potential Vij{r) obtained with a = oo becomes identical to the potential 
Vij{r) corresponding to the case a = 0. Hence the the case a = oo appears to be identical to the 
case a = at any energy E including E = Eq. To demonstrate this explicitly, let us examine the 
wave function Lpi{EQ,r) in the case a = oo. Substituting E by Eq in (|18D we obtain: 



ifi{Eo,r) = (21) 

j oo 



Aipi{Eo,r) 



A-C + CZl\vAEo,s)\' ds 
j 



(22) 



Replacing Ahj A + C and normalizing the wave function (^2]), we obtain the expression (|T 

B. Supersymmetry 

Let us discuss a particular choice of parameters: C = 1, a = oo, and A = 1. The wave function in 
this case is 

oo 

ipi{Eo,r)^J(Pj*{Eo,s) ^j{E,s) ds 
^E, r) = ip.iE, r) + '-^ ^ (23) 

j 

r 

ViiEo, r) X; / Vj*iEo, s) ^j{E, s) ds 
= VriE, r) 1^ ^ . (24) 

j 

Equation (^3]) can be used at any energy E while Eq. (|2^) is applicable only ii E ^ Eq. In the case 
E = Eq, the wave function can be rewritten in a simpler form as 

(Pi{Eo, r) = . (25) 

'Zl\'PjiEo,s)f ds 
j 



Equation ( [23| ) is just the Eq. (4) of Ref. |2y]. In Ref. |]20|, Sparenberg and Baye suggested a 
multi-channel supersymmetry transformation. Thus equations (p^)-(|25|) describe the multi-channel 
supersymmetry transformation, or, in other words, the multi-channel supersymmetry transformation 
is a particular case of the phase-equivalent multi-channel transformation discussed in this paper that 
corresponds to the particular choice of the parameters. Let us discuss how it works. 

It is clear from Eq. ( p5|) that \ipi{Eo, r)| — ^ oo as r ^ 0. Hence at the energy Eq, the wave function 
ipi{EQ,r) does not match the required boundary condition (^) at r = 0. At the same time, ipi{Eo,r) 
fits the boundary condition (^ at r = oo. Therefore it is impossible to construct another solution 
of the Schrodinger equation (0) consistent with both boundary conditions at the energy E = Eq. 
As a result the phase-equivalent transformation removes the bound state at E = Eq. At the same 
time, it is clear from (p^ that for all energies E ^ Eq, zero in the denominator arising in the limit 
r — >■ is canceled by the zero in the numerator and the wave function (^) matches the boundary 
conditions at the origin and at the infinity both at once. So, the transformation in this case removes 
the bound state at E = Eq but not any of the other bound states, while the S'-matrix at any energy 
E > is unchanged. 



Of course, the supersymmetry transformation can be also formulated in the case a = 0. It is 
interesting that the bound state in this case is removed by a different mechanism. Suppose A = 
and C is arbitrary. The wave function at any energy E in this case may be written as 

r 

^i{Eo,r)^JiPj*{Eo,s) ^j{E,s) ds 
if^iE, r) = ^,{E, r) '-^ ^ . (26) 

3 

It is seen that aX E = Eq the wave function (fi{EQ, r) = 0. 

We used the boundary condition (^ to construct the supersymmetry transformation: the bound 
state is removed because for some particular parameter values the wave function (^j(i?o, r) diverges at 
the origin and appears to be inconsistent with (^). One can suppose that it is also possible to use the 
boundary condition at r = oo to remove the bound state and to construct another supersymmetry 
transformation. It is not so. Let us discuss the case oi a = oo, A = and arbitrary C. As is seen 
from (P^), ipi{EQ,r) = in this case, thus the bound state is removed. However the transformation 



is no more phase-equivalent. Really, at energies E > Eq the last term in (|T8|) does not vanish at 
r — > cxD and provides an additional phase shift, or, in other words, it modifies the S'-matrix. 
An example of an application of the multi-channel supersymmetry transformation to the Moscow 



A^A^ potential can be found in Ref. ||20|| . 



C. Inverse supersymmetry 

A transformation that adds a bound state to the discrete spectrum of the system and leaves 
unchanged the S'-matrix and the energies of all bound states supported by the initial Hamiltonian, 
we shall refer to as inverse supersymmetry transformation. 

Let us suppose that there is no bound state at the energy £"0 < 0. By Lpi{EQ,r) we now denote 
the wave function at energy Eq that matches the boundary condition (0) at infinity but diverges at 
the origin as r~'' (see, e.g.Q, |]2I|) where /« is the angular momentum in the channel i. 

Having ipi{EQ,r), we can use our transformation to obtain the homogeneous Schrodinger equation 
(p^) in the case a = 00. The transformed wave function (fi{E,r) is given by (|15). It is seen from 
(|T^) that Lpi{E,r) does not diverge in the origin and matches the boundary conditions both at 
the origin and at infinity, at any energy E ^ Eq. For E = Eq, the transformed wave function 
(pi{EQ,r) is given by (^l]) . It is clear that (pi{EQ,r) at the origin is proportional to r^^~''~^ where 
L = max{/j}. Hence, (pi{EQ,r) matches the boundary condition (^ if L > 2 and is not consistent 
with (^) if L < 1. Therefore our transformation with ifi{EQ, r) irregular at the origin, is the inverse 
supersymmetry transformation in the case L > 2. In the case L < 1 the transformation appears 
to be a phase-equivalent transformation that does not make use of the bound state and can be 
applied to the system that does not support a bound state. If the transformation is applied to the 
free Hamiltonian with Vij{r) = in the s or p partial wave, it produces a non-zero 'transparent' 
potential Vij{r) that provides phase shift 5 = at any energy E. The multi-channel version of the 



'^The r~ divergence of the wave functions at the origin is derived in Ref. pl| for the single-channel case only. However the 
derivation of the r~ rule of Ref. ^Ij] can be easily generalized to the multi-channel case, at least for the potentials that do not 
diverge in the origin. 



transformation couples s and p partial waves to produce a two-channel 'transparent' interaction that 
provides the S'-matrix of the form Sij = 6ij . 

It is interesting that the inverse supersymmetry transformation is not unique: we have three 
parameters Eq, A and C that provide a family of inverse supersymmetry partner potentials. Contrary 
to it, the supersymmetry transformation is unique; however, it can be used in combination with the 
phase-equivalent transformation to construct a family of potentials phase-equivalent to the initial 
one but not supporting one of the bound state. 

The multi-channel inverse supersymmetry transformation is discussed in more detail in a very 
recent paper of Leeb et al |2^. One can find in this paper examples of the applications of the 
transformation to realistic A^A^ potentials. This transformation is discussed in Ref. |^, too; in 
particular the authors of Ref. [^ also conclude that it is possible to create a new bound state by 



means of the phase-equivalent transformation only in the case L > 2. 

IV. CONCLUSIONS 

We derived a multi-channel phase-equivalent transformation that can be used without restric- 
tions on the structure of the discrete spectrum of the system in various scattering problems like 
A^A^ scattering, nucleon-cluster or cluster-cluster scattering. The multi-channel supersymmetry and 
inverse supersymmetry transformations appear to be particular cases of the suggested general phase- 
equivalent transformation corresponding to particular choices of the parameter values. The inverse 
supersymmetry transformation is possible if only the orbital angular momentum /» > 2 at least 
in one of the coupled channels. It is interesting to note that from the point of view of the A^A^ 
system, this means that a deep attractive A^A^ potential supporting an additional forbidden state 
like Moscow A^A^ potential, can be constructed by the inverse supersymmetry transformation of 
the realistic meson-exchange potential with repulsive core only due to the d wave admixture in the 
deuteron wave function. 

With the help of the suggested transformation, one can construct a family of phase-equivalent 
potentials depending on continuous parameters. Such families may be very useful for fine tuning of 
the interaction aimed to fit not only two-body observables but also three- and few-body ones. If 
the system has at least one bound state, the phase-equivalent potential family is constructed using 
directly formulas (H) and @. One can construct phase-equivalent single- or multi-channel potential 
families also in the case when there are no bound states in the system: if all channel orbital angular 
momenta li < 1, one can apply directly the transformation with the irregular function ipi{EQ,r)] if 
at least one of the channel orbital angular momenta li > 2, one can produce a bound state using 
inverse supersymmetry at the first stage and remove the bound state at the last stage with the help 
of supersymmetry version of the transformation. So, one can, for example, construct a family of 
phase-equivalent potentials for any combination of coupled partial waves in the A^A^ system. 

We hope that the suggested transformation will be useful in various few-body applications. 
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